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Abstract 

The inverse transfer in two-dimensional turbulence governed by the surface quasi- 
geostrophic (SQG) equation is studied. The nonlinear transfer of this system con- 
serves the two quadratic quantities ^'i = (|(-A)^/'^'i/'P)/2 and ^'2 = (|(-A)^/2^/;p)/2 
(kinetic energy), where ip is the streamfunction and (•) denotes a spatial average. 
In the limit of infinite domain, the kinetic energy density ^'2 remains bounded. For 
power-law inverse-transfer region, the inverse flux of ^'i diminishes as it proceeds 
toward sufficiently low wavenumbers, implying that no persistent inverse cascade of 
^1 is sustainable. The unrealizability of an inverse cascade of ^'i implies that there 
is no direct cascade of ^'2- Hence, the dual-cascade picture which is widely believed 
to be realizable in two-dimensional Navier-Stokes turbulence does not apply to SQG 
turbulence. Numerical results supporting the theoretical predictions are presented. 
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1 Introduction 



The motion of a three-dimensional stratified rotating fluid is characterized 
by the geostrophic balance between the Coriolis force and pressure gradient. 
The dynamics governed by the first order departure from this linear balance 
is known as quasi-geostrophic dynamics (see for example [4,5,20,22]), which 
can be described in terms of the (three-dimensional) geostrophic streamfunc- 
tion ^/'(x, t). The vertical dimension z is usually taken to be semi- infinite and 
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doubly periodic conditions are usually imposed on the horizontal flow. Nor- 
mally decay conditions are required as z ^ oo. At the flat surface boundary 
z — 0, the vertical gradient of ip{'x,t) matches the temperature field T(x, t), 
i.e. T(x, t)|^=o = (^^^(x, t)|^=o- For flows with zero potential vorticity, this 
surface temperature field can be identified with {—Ay^^t/j, where A is the 
(horizontal) two-dimensional Laplacian. The conservation equation governing 
the advcction of the active temperature field (—Ay^'^ip by the surface fiow 
u={-dyilj,d,ilj) is [2,13,20,21] 

a,(-A)V2^ + J(V',(-A)V2^) = 0, (1) 

where J{-d, 9) = d^-^dyO — d^Ody-d. Eq. (1) is known as the SQG equation. 

In this paper, a forced-dissipativc version of (1) is studied. A dissipative term 
of the form fiAip, where /i > 0, which results from Ekman pumping at the 
surface, is considered (cf. [6,24]). Since (— A)^/^-?/^ is the advected quantity, 
this physical dissipation mechanism corresponds to the hypoviscous dissipation 
operator A)^/^. The dissipation coefficient jj, has the dimensions of velocity 
and is not vanishingly small in the atmospheric context. The system is assumed 
to be driven by a forcing /. Thus, the forced-dissipative SQG equation can be 
written as 

a* (-A)^/V + Ai^, (-A)^/V) = /^Av^ + /■ (2) 

The Jacobian operator J(-, •) admits the identities 

(0J(^, 6)) = -(^J(0, 6)) = -{em 0)), (3) 

where (■) denotes the spatial average. As a consequence, the nonlinear term 
in (2) obeys the conservation laws 

(^J(^, (-A)VV)) = ((-A)VVJ(^, (-A)VV)) = 0. (4) 

It follows that the two quadratic quantities = (|(— A)-'^/^'0p)/2 and ^2 = 
(|(— A)^/^'0p)/2 (kinetic energy) are conserved by nonlinear transfer. 

The simultaneous conservation of ^'i and ^2 by advective nonlinearities im- 
poses strict constraints on the transfer (redistribution) of these quantities in 

wavenumber space. Let us consider the transfer of an amount ^1 = e, initially 
distributed around a given wavenumber s, which corresponds to an initial 
^2 = se. Let 'ifi{k) and ^'2(/c) be the resulting redistributions of ^'i = e and 
of ^^2 — se. Given arbitrary wavenumbers p < s and q > s, one has 
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- /*i(A;)dA;<— /*2(A;)dA; < -, (5) 
e J qe J q 

Q g 

p p 

— I ^!2{k)dk<— I ^!Ak)dk<-, (6) 
se J se J s 



where the two conservation laws and straightforward inequahties have been 
used. The left-hand side of (5) [(6)] is the fraction of e [se] that gets trans- 
ferred to wavenumbers k > q [k < p\. This fraction is bounded from above 
by s/q [p/s], implying that no significant fraction of e [se] can be transferred 
to wavenumbers A; ^ s [A; <S s]. This type of constraint on the nonlinear 
transfer of the invariants is a common feature in incompressible fiuid systems 
in two dimensions. (Some familiar systems in this category are the Charney- 
Hasegawa-Mima equation [11,12] and the class of a turbulence equations [21], 
which includes both the Navier-Stokes and the SQG equations.) The implica- 
tion is that when the said initial sources spread out in wavenumber space, 
[^2] is preferentially transferred toward lower [higher] wavenumbers. Accord- 
ing to the classical theory of two-dimensional turbulence [1,14,15,17], which 
was originally formulated for high-Reynolds number Navier-Stokes fiuids and 
subsequently thought to apply to two-dimensional incompressible fluids in 
general, this preferential transfer achieves the extreme limit by transferring 
virtually all e to /c <^ s (inverse cascade) and virtually all se to /c ^ s (direct 
cascade). The transfer of the invariants in this manner is known as the dual 
cascade. ^ 

However, it is shown in this work that the preferential transfer of ^1 and of 
^^2 is not as dramatic as predicted but rather has a limited extent. More ac- 
curately, it is shown that an upper bound on the inverse flux of across a 
low wavenumber ^ vanishes (uniformly in time) as £/s — » 0, where s is the 
characteristic forcing wavenumber, thereby ruling out the existence of a per- 
sistent inverse cascade of ^1 toward the low wavenumbers. The unrealizability 
of an inverse cascade implies that there is no direct cascade of \E'2 [24]. The 
physical reasons behind this behaviour are that the kinetic energy density of 
SQG dynamics remains bounded in the limit of inflnite domain, a consequence 
of the hypoviscous dissipation operator /i(— A)^/^, and that the inverse flux of 
^1 across I is proportional to the energy content of the wavenumber region 
k < i. These two properties of SQG turbulence, when combined, imply that 
the inverse flux of ^1 becomes smaller for progressively lower £, thus ruling 
out the existence of a persistent inverse cascade. 

In the next section some preliminary estimates, which are employed in the 



For some recent discussion on the possibility of a dual cascade in various 
two-dimensional systems, including the Navier-Stokes and SQG equations, see 
[24,25,27,29]. 
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derivations of the main results in Section 3, are presented. Section 4 examines 
the imphcations of the results in Section 3 for the long-time dynamics and 
spectral distribution of kinetic energy. Section 5 presents some numerical re- 
sults in support of the theoretical prediction of no cascades in SQG dynamics. 
The paper ends with some discussion in the final section. 



2 Preliminciry estimates 



This section presents a simple analytic inequality and reviews the boundedness 
of the energy density [24] . The former applies to bounded domains only, while 
the latter is valid for both bounded and unbounded cases. 

For a doubly periodic domain L x L, the (complex) Fourier representation of 

ip is 



V'(x)=^exp{zk-x}^(k). (7) 

k 

Here k = kQ{n,m), where = 2ti / L is the lowest wavenumber and n and m 
integers not simultaneously zero. For a given wavenumber let ip'^ and ip^ 
denote, respectively, the components of ip spectrally supported by the disk 
d — {\^\k<^} and its complement D = {]i\k> ^}, i.e. 

-0^ = ^ exp{ik • x}'^(k), V''^ = X/ 6xp{ik • x}'^(k). (8) 

ked kGD 

For the lower- wavenumber component V'^, the following inequality holds: 



1/2 / 1/2 

,<V2 



sup|vn<E^wk)i< Em \T.k'\m?\ -c-^f' (9) 

kGd VkGd / VkGd / 

where c is an absolute constant of order unity and ^2 — (I^V'^P)/2 the 
large-scale energy density associated with the wavenumbers k < (9) the 
Cauchy-Schwarz inequality is used in the second step, and the sum X^ked 1 ~ 
/kl represents the number of wavevectors in d. 

On multiplying (2) by ip and (— A)^/^-;/' and taking the spatial averages of the 
resulting equations over the domain, noting from the conservation laws that 
the nonlinear terms identically vanish, one obtains evolution equations for 
and ^2, 
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^*i = -2/iVl/2 + (M, (10) 
^*2 = -2/.*3 + (/(-A)VV), (11) 

where ^3 = (|(— A)^/^'0p)/2. Using the Cauchy-Schwarz and Young inequal- 
ities, one obtains the upper bounds on the injection terms in (10) and (11): 



(/(-A)V2^) < (|(-A)-VVp)i/^(|(-A)-VVHi/^ < + (12) 

where the norms of / are defined by F_i = (|(— A)^^/'^/p)/2 and F_2 = 
(|(— A)"^/^/P)/2. These norms are assumed to be finite. Substituting (12) in 
(10) and (11) yields 



< -i^-^2 + l^-'F_2, (13) 

-^*2<-//*3 + /^"'F-i. (14) 
dr 

Let the overline denote the asymptotic average (existence is assumed). One 
can deduce the following upper bounds for ^2 and ^3: 



^2</X-'^-2, (15) 

%<l^-^F-i. (16) 



Ineqs. (15) and (16) apply to both the unbounded and bounded domains. In 
the study of turbulent transfer, the spectral support of / is required to be 
in the intermediate wavenumber region, so as to render inverse- and direct- 
transfer ranges free of sources. The hmit of infinite domain {ko — > 0), for which 
the classical theory is formulated, is taken in a straightforward manner: in the 
limit L 00, the dissipation coefficient n, the injection densities, and the 
forcing characteristic length scale (cf. [30]) are held fixed. In this limit both 
F_i and F_2 are bounded if F = (|/p)/2 is bounded since F_2 < k^F_i < 
k^F, where k^ is the minimum wavenumber of the spectral support of /. A 
persistent inverse cascade of toward ever- lower wavenumbers necessarily 
requires d^i/dt > 0, which, by (13), implies ^2 < ^^"^-^-2- Therefore, since 
the main concern of the subsequent analyses is the realizability of an inverse 
cascade, for the rest of this paper the (instantaneous) energy ^'2 is assumed 
to be bounded in the limit ko ^ 0. 
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3 Diminishing inverse transfer 



This section reports the main result of this paper. In the subsequent analyses, 
rigorous estimates are supplemented by the usual assumption of power-law 
scaling for the energy spectrum. An inverse-transfer range '^2{k) = ak'"', for 
A; < s, is assumed. Here s may be taken to be the minimum wavenumber of the 
spectral support of the forcing. An upper bound for the total transfer of \l'i 
into the low- wavenumber region [ko,i <S s], i.e. the inverse flux of ^'i across 
£, is derived in terms of the kinetic energy ^2 and of the large-scale kinetic 
energy It is shown that this flux diminishes in the limit i/s — >■ 0. The 
implication is that no persistent inverse cascade of ^1 is realizable. 

The evolution of *f = (|(-A)V4^<|2^/2 is governed by 



= -(V'<J(V>,(-A)VV))-2/.*< 
= ((-A)V2^j(V,>,V'<))-2/.*< 

<(|(-A)VV||VV'>||Vn) -2/.*2<, (17) 

where (3) has been used in both the second and third steps. The final triple- 
product term represents an upper bound for the inverse flux of ^1 across the 
wavenumber i that drives the large-scale dynamics. 

In the limit of infinite domain, no finite rigorous estimates of the nonlinear 
term in (17) are available. Nevertheless, since the norm [V^'^l represents a 
measure of the large-scale fluid velocity (associated with k < £), it could 
heuristically be identifled with ^2^^^. Hence, a rough estimate of |V'0^| would 
be |V^<| ^ 2-i/2c'(|V^<p)^/2 = c'*|^^/^ where c' is a constant. Substituting 
this estimate into (17) yields 



-^*< < 2c'*<'/'*2/'*>'/' - 2//*< < 2c'*<'/'*2 - 2/x*<, (18) 
dt 

where = {\'Vip^\'^) /2 is the small-scale energy, associated with k > £. 
Since \E'2 is bounded and since its spectrum '^2{k) is assumed to obey power- 
law scaling in the inverse-transfer region, ^'^ necessarily diminishes as £ — > 0. 
Therefore, the quantity 2c'^2 ^''^\E'2, which bounds the inverse flux of ^1 across 
i, can become arbitrarily small for sufficiently low i. Thus, no persistent inverse 
cascade of \&i is realizable provided that the foregoing heuristic estimate of 
IVV'^I can be assumed. 

A rigorous version of the above calculation can be deduced for the bounded 
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case. By applying (9) to (17) one obtains 



;^vi/<<sup|V^<|(|(-A)^/V||VV'>|)-2/.^< (19) 
at X 

< 2cf *f /^^f *>^/^ - 2;.^< < 2cf (^] - 2;.^<. 

The difference between (18) and (19) is tlie presence of the ratio i/ko in the 
latter. Hence, the arguments in the preceding paragraph go through without 

change if the ratio i/ko is either held fixed or allowed to grow at some certain 
rate as £ — * 0. This condition can be stated more explicitly by making use 
of the assumed spectrum ^'2(/c) = aA;~° for < s. In the limit ko — > 0, the 
large-scale energy and the total energy ^2 can be estimated as 



= a y A;-" dA; = j^^'"", (20) 
^ 

s 

■^2>a f k-"dk= -^s^"^, (21) 
J 1 — a 








where a < 1, in accord with the boundedness of energy. It follows that 



*f<2cl(g »f. (22) 

Given a bounded ^2 and a ratio l/ko that does not grow as rapidly as 
(s/iy^~"^^'^ as £ becomes small, this upper bound on the inverse flux of ^E^i 
clearly diminishes as i/s ^ 0. Hence, for a sufficiently wide inverse-transfer 
range [i,s], the advective nonlinearities of SQG turbulence are incapable of 
transferring a significant amount of the injection of ^'i to the low-wavenumber 
region [ko,£]. 

It is interesting to generalize the above calculations to other models of incom- 
pressible fluid turbulence in two dimensions. Pierrehumbert et al. [21] propose 
to consider the so-called a-turbulence models, for which the unforced and 
inviscid dynamics are governed by 



9t(- A)°/2V + A^, (-A)"/^) = 0, (23) 

where « is a positive number. The two invariants of this system are \E'q = 
(|(_A)°/4^|2)/2 and ^2a = (|(-A)"/VP)/2. The inverse transfer of across 
i can be estimated as 
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I (V'< J(V', (-A)"/^V)) I = I ((- A)-/^^' I 

<(|(-A)'^/V||VV'>||Vn) 
<sup|VV'<|(|(-A)"/V||VV'>|) 

X 

<24*f/^<^*>^/l (24) 

There are cases for which ^'^ — > in the hmit £ — > 0, leading to the di- 
minishing of the inverse flux across i. For example, if the direct-cascading 
candidate \E'2q remains bounded^ in the limit ko ^ and if a < 1, then 
^< < £2-2a^<^ ^ £2-2a^|(_^)a/2^<|2^/2, which Certainly vanishes as f 
because ^2a) is bounded. On the other hand, ^2 converges toward the 

low wavenumbers for the same reason. There remains the modest assumption 
that ^2 ^Iso converges toward the high wavenumbers. Given all these, the 
preceding arguments of no inverse cascade go through without change. 

Remark 1. For technical reasons, it is difficult to generalize the analyses in 
this section to the truly unbounded case (i.e. the case £/ko —>■ 00, regardless of 
the ratio i/s). The main difficulty lies in the intrinsic domain-size dependence 
of the nonlinear term, thereby making its rigorous estimates, such as the one 
in (19) and the subsequent estimate (22), diverge in the limit kg 0. 

Remark 2. If the viscous operator /i(— A)^/^ is replaced by one of higher 
degree, as is often done in numerical simulations, then the boundedness of 
energy in the limit ko ^ cannot be guaranteed. Nevertheless, if boundedness 
of energy is assumed, then the above calculations go through without change 
since they do not refer at all to the dissipation mechanism. 

Remark 3. The above procedure, when applied to the nonlinear term of the 
two-dimensional Navier-Stokes equations, yields an estimate of the inverse 
energy flux that is energy-dependent. This behaviour makes it challenging to 
estimate the inverse energy transfer since the energy is supposed to grow with 
progressively wider inverse-transfer range. Nevertheless, it can be shown that 
for the Kolmogorov-Kraichnan k~^^^ energy spectrum the energy that gets 
transferred onto ko can be bounded from above by a constant independent of 
the width of the inverse-cascading range, i.e. independent of A'o. This result, to- 
gether with related issues such as the Kolmogorov constant and the Kraichnan 
conjecture of energy condensation at ko, is the subject of a separate study. 

Remark 4. Recently, Eyink [10] raised the possibility of energy condensa- 
tion in bounded two-dimensional Navier-Stokes turbulence, a conjecture by 
Kraichnan [14], and argued that this possibility cannot be ruled out in some 

^ This is certainly the case for SQG (a = 1) and Navier-Stokes (a = 2) turbulence, 
for which the energy and enstrophy, respectively, remain bounded in the limit of 
infinite domain. 
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recent theoretical results. Kraiclinan [14], in an attempt to apply his dual- 
cascade hypothesis (initially formulated for unbounded fluids) to bounded tur- 
bulence, predicts that the inverse energy cascade, upon reaching ko, deposits 
energy to this wavenumber, and that this process continues until growth of 
energy at ko is limited by its own dissipation, resulting in what may be termed 
an "energy condensate". The energy dissipation by this condensate alone is 
supposed to account for virtually all the energy injection, so that the energy 
condensate is also an enstrophy condensate, although the latter is of a lesser 
degree. The realization of such a "singular" energy and enstrophy concentra- 
tion at ko (or around ko) is required to maintain the proposed dual cascade 
in the bounded case. Recent numerical results seemed to suggest otherwise: 
as the turbulence approaches a steady state, a k^^ energy spectrum forms at 
the large scales [3,27]. Nevertheless, the a priori exclusion of the Kraichnan 
scenario by some recent theoretical studies, such as Constantin et al. [7], Tran 
and Shepherd [29], and Kuksin [16], does not seem to be fully justified. The 
SQG dynamics allows for no possibihty of such a condensate. 



4 Approach to steady dynamics and spectral distribution of energy 

This section features some physical interpretations of the results derived in the 
preceding section. The non-cascading dynamics of SQG turbulence is discussed 
together with a review of the constraint on the spectral distribution of energy 
derived by Tran [24]. 

It is customary in the study of 2D turbulence to consider the scenario in which 
the fluid is driven around a wavenumber s by steady injections (fip) = e and 
(/(— A)"'^/^-?/') = se. The result in the preceding section implies that an inverse 
transfer of a nonzero fraction of e to sufficiently low wavenumbers requires that 
the energy spectrum in the inverse-transfer region be no shallower than ak~^. 
But then the energy would grow at least as rapidly as aln(s/£) as s/i ^ oo, 
eventually leading to a balance between the injection e and the dissipation 
2/i\l'2. This result suggests two plausible routes to the long-time high- Reynolds 
number dynamics (for some suitably defined Reynolds number). First, if a k~^ 
inverse-cascading range is realized, ^ the inverse flux decreases as it proceeds 
toward lower wavenumbers since the dissipation of ^'i, given by 2/x^2, grows 
logarithmically, eventually becomes steady, as described above. Second, 
suppose that an inverse-cascading range with energy spectrum shallower than 
k~^ is realized. The inverse flux decreases for the reason discussed in the 
preceding section as it proceeds toward sufficiently low wavenumbers. As a 
result, growth of (and of ^^2) occurs alongside the existing inverse-transfer 
range. Eventually the dissipation 2/i^2 reaches the injection e and ^1 becomes 

^ Incidentally, dimensional analyses predict a inverse-cascading range. 
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steady. If power-law scaling is maintained, the slope of the energy spectrum 
in the inverse-transfer range approaches — 1 . The low- wavenumber end of this 
range £ can be calculated from 2/ialn(s/£) = e. 

Near steady dynamics can be achieved after the inverse flux of across 
I becomes sufficiently less than its dissipation 2yU^2) which should then be 
comparable to e. Hence, by replacing \E'2 in the expression on the right-hand 
side of (22) by e/2/x and requiring that the resulting expression be no larger 
than e, one obtains a condition for this near steady picture: 




where a constant factor of order unity has been dropped. For fixed (per-unit- 
area) injection rate e, forcing wavenumber s, and dissipation coefficient ^, no 
significant fraction of e can be transferred to wavenumbers k < £, where £ 
satisfies (25). ^'i necessarily becomes near steady, with the low- wavenumber 
region [/cq, i] at best weakly excited. 

For steady dynamics the balances 2/i\l'2 = e and 2/i^'3 = se are achieved. It 
follows that s^2 = ^3, or in terms of the energy spectrum ^'2(^), 

oo 

J (s - k)'^2(k) dk ^ 0. (26) 

feo 

This equation can be used to estimate the slopes of the energy spectrum if 
power-law scaling is assumed for both the inverse- and direct-transfer ranges. 
For this purpose, let us consider the following spectrum 

"i<l<l, 

where a, b, a, (3 are constants, and kij is the highest wavenumber in the 
range k^l^ , beyond which the spectrum is supposed to be steeper than k^^ . By 
substituting this spectrum into (26) and making the respective substitutions 
K = k/s ioT k < s, and k = s/k for k > s, one obtains [24] 

1 1 

j {1- KjK-'^dKP^ j {l-K)K^-^dK, (28) 

t/s s/kt, 

where the contribution from both k < i and k > k^ has been dropped. It 
follows that \l l/s > s/kv, then —a < /3 — 3. Hence, the constraint 
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a + (3>3 



(29) 



holds. Since a < I, P satisfies /? > 2, meaning that for k > s, the spectrum 
^3(A;) of \l/3 is no shallower than k~^. Hence, the energy dissipation cannot 
occur mainly at A; ^ s. Thus there is no direct cascade, a dynamical behaviour 
consistent with no inverse cascade of ^'i. 

Remark 5. Although no persistent inverse flux is possible, in the limit ko — > 0, 
^1 could become unbounded, growing at a rate that fluctuates about zero and 
has a vanishingly small positive average. The unboundedness of ^'i requires 

that the energy spectrum of the inverse-transfer range have a non-positive 
slope, so that the spectrum '^i{k) is at least as steep as kr^. For a given set 
of physical parameters, the issue of whether or not \l'i becomes divergent (in 
the limits /cq — > and t ^ oo) is interesting but is beyond the scope of the 
present work. 

Remark 6. In some sense SQG turbulence is relatively "simple" . Given spec- 
trally localized steady injections about the forcing wavenumber s, the dy- 
namics should eventually become non-cascading. For a k~^ transient inverse- 
cascading range, the approach to steady dynamics is rather slow for two rea- 
sons. First, because ^i(A';) oc A;^^. the low- wavenumber end of the inverse- 
cascading range i proceeds relatively slowly toward lower wavenumbers, even 
during the early stages for which the inverse cascade is relatively strong 
{d'^i/dt e). Second, growth of energy is only logarithmic in giving rise 
to rather slow growth of the dissipation of toward the lower wavenumbers. 

Remark 7. Dimensional analyses, without references to any particular dis- 
sipation mechanisms, predict a k~^ inverse-transfer range and a k^^^^ direct- 
transfer range. The former is consistent with an inverse cascade of ^i, but the 
latter allows for virtually no energy to get transferred to the small scales. ^ 
Before a viscous dissipation mechanism is taken into consideration, a k^^^^ 
spectrum means that virtually no energy gets transferred to the high wavenum- 
bers. In the presence of a viscous dissipation operator of the form oc (—A)'', 
for < 5 < 1/3, instead of the natural dissipation operator //(— A)^/^, a 
^-5/3 (direct-transfer range means that the spectral energy dissipation scales 
as A;"^/^"*"^^, which is no shallower than k~^, thereby allowing for virtually no 
energy to be dissipated at its high-wavenumber end. 

Remeirk 8. The simultaneous conservation of and leads to an increase 



These predictions seem to be mutually inconsistent since the inverse-transfer of 
*i via a spectrum is incompatible with the "frozen-in" of energy due to the 
^-5/3 direct-transfer range. This picture is quite contrary to that of two-dimensional 
Navicr-Stokes turbulence, for which a k^"^ enstrophy-transfer range means that 
virtually all enstrophy gets transferred away from the forcing region, even before 
dissipative effects are considered. 
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in ^4 (enstrophy) when an initial spectral peak spreads out in wavenumber 
space [24] . This explains the observed formation of strong "fronts" in numerical 
simulations of SQG turbulence [8,9,18], even in the absence of a direct energy 
cascade as discussed above. 



5 Numerical results 



This section reports results from numerical simulations that illustrate the 
diminishing inverse transfer and no cascades of SQG dynamics. Numerical 
studies in the literature have thus far failed to recognize these properties of 
SQG turbulence (see for example [19,21,23,26,28]). Tran and Bowman [28], 
however, notice that is "reluctant" to cascade to the large scales, even 
when the natural dissipation operator A)^/^ is replaced by ones with higher 
degrees, allowing for relatively weaker dissipation at the large scales. 

Equation (2) is simulated in a doubly periodic square of side 27r, where the 
modal forcing /(k) is nonzero only for those wavevectors k having magnitudes 
lying in the interval K = [9.5, 10.5]: 



/(k) = ^ ^ , (30) 

^^EipmIV'Ip)!^ 

where se = 1 is the constant energy injection rate and the number of 
discrete wavenumbers in K. The wavenumber s ~ 10 is defined such that s"^ 
is the mean of over K. The (constant) injection rate of is e ~ 0.1. 
Dcaliascd 683^ and 1365^ pseudospcctral simulations (1024^ and 2048^ total 
modes) were performed. For Navier-Stokes turbulence, these resolutions are 
sufficient to simulate an inverse cascade that carries about a quarter of the 
energy injection to the large scales via a discernible k~^^^ range [24]. For SQG 
turbulence, even the higher resolution turns out to be insufficient to simulate 
a noticeable transient inverse cascade. Two dissipation coefficients were used: 
/J, = 0.05 and /i = 0.025. The lower and higher resolutions were used for the 
stronger and weaker dissipation, respectively. Both simulations were initialized 
with the spectrum *2(^) = W'^nk / {100 + k^). 

Figure 1 shows the time-averaged (between t = 19.7 and t = 20.3) ^ kinetic 
energy spectrum for the case //. = 0.05. The average energy is ^'2 = 0.99. 
The dissipation of ^1, averaged for the same period, is 2ix^2 = 0.099, which 



^ The dissipation time at the forcing wavenumber s ~ 10 is {2fis)^^ ~ 1, so that 
the time t ~ 20 is sufficiently long for the evolution of ^'2. In fact, it was observed 
that both *i and ^^2 became steady at t pa 12. 
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Fig. 1. The time-averaged steady-state energy spectrum ^'2(A:) vs. k for the dissi- 
pation coefficient fi = 0.05. The average energy is *2 = 0.99, implying that the 
dissipation of ^'i, averaged in the same period, is 0.099. This amounts to virtuahy 
ah of the injection rate e 0.1. Hence, no inverse cascade of ^'i exists and both 
and ^'2 are steady. 

amounts to virtually all the injection rate ~ 0.1. The energy was observed 
to increase monotonically from t = to t ^ 12, hj which time d\l/i/dt ~ 0; 
no significant inverse transfer was observed throughout the period. Instead, 
growth of energy occurs mainly around the forcing region. The steady spec- 
trum is relatively shallow in the low-wavenumber range (see Fig. 1), meaning 
that the lowest wavenumbers are virtually unexcited. The small-scale energy 
spectrum scales as k^^'^, so that the spectrum ^3 (A;) of the energy dissipation 
agent ^'3 scales as A;~^'®. This scaling means that the energy dissipation occurs 
mainly around the forcing region, consistent with the weak inverse transfer. 

A somewhat stronger transient inverse transfer was observed for the case = 
0.025. Fig. 2 shows a near steady kinetic energy spectrum averaged between 
t = 15.1 and t = 15.6. The average energy is ^2 = 1-96. The dissipation of 
^1 averaged for the same period is 2/x\I/2 = 0.098, which amounts to virtually 
all of the injection rate f« 0.1. The large-scale spectrum is better "filled up" 
than that of the previous case, due to a stronger inverse transfer during the 
transient phase. Nevertheless, no significant fraction of e reaches the lowest 
wavenumbers. The energy dissipation occurs around the forcing region, as is 
evident from the steep small-scale spectrum (see Fig. 2). 

It is expected that for higher resolutions (so that simulations with smaller fi 
are possible), transient inverse fiuxes can become more noticeable and steeper 
inverse-transfer ranges can be realized. Given limited resolutions, Tran and 
Bowman [28] use high-degree dissipation operators oc A and oc (— A)^/^. In 
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Fig. 2. The time-averaged near steady-state energy spectrum '^2ik) vs. k for the 
dissipation coefficient n = 0.025. The average energy is ^'2 = 1.96, implying that 
the dissipation of averaged in the same period, is 0.098. This amounts to most 
of the injection rate e ~ 0.1. 

both cases, weak inverse cascades are observed, but the transient energy spec- 
tra of the inverse-transfer region are considerably shallower than k~^. These 
spectra cannot support a persistent inverse cascade, as argued in Section 3. 
However, an interesting possibility arises. Due to both the inability of the 
nonlinear transfer to sustain a constant (wavenumber-independent) inverse 
flux and d^i/dt > 0, growth of ^1 ought to occur within the existing inverse- 
transfer range, thereby causing this range to become steeper. This process may 
continue until the inverse-transfer range eventually exceeds the k~^ thresh- 
old. Because of the high degrees of viscosity, such a slope steepening causes 
insignificant increase in the dissipation of ^1. As a result, a positive growth 
rate d^i/dt could be maintained for a steeper-than-A;~^ inverse-transfer range, 
which might then be able to support a persistent inverse cascade. 



6 Conclusion 

In this paper, the advective transfer of SQG turbulence is studied. The main 
result obtained is an upper bound for the inverse flux of the inverse-cascading 
candidate ^1. This upper bound diminishes as the flux proceeds toward suf- 
ficiently low wavenumbers, thereby ruling out the existence of a persistent 
inverse cascade. The unrealizability of an inverse cascade entails that there 
is no direct cascade. This is the first rigorous example of the dynamics of 
incompressible fluids in two dimensions that exhibits no cascades. 
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There are two essential features of SQG turbulence that facilitate the proof of 
non-cascading dynamics. First, the inverse flux of across a low wavenumber 
£ can be uniformly (in time) estimated in terms of the kinetic energy \I'2 and of 
^2 ) the large-scale energy associated with the low- wavenumber region k < i. 
Second, the energy ^2 remains bounded in the limit of infinite domain. The 
former is an intrinsic property of the advectivc nonlinear term, and the latter 
is due to the hypoviscous nature of the dissipation of SQG dynamics. The 
boundcdness of energy means that if the turbulence is driven at some fixed 
energy density rate around some fixed wavenumber s, virtually no energy can 
be acquired by wavenumbers k <^ s. This means that in the limit i/s — > 0, 
^2 can become arbitrarily small, leading to an arbitrarily small inverse fiux 
of \&i across £. Hence, for sufficiently low wavenumber £ no significant inverse 
transfer of ^'i across i is possible. 

Numerical simulations of SQG turbulence were performed. The results show 
no significant inverse transfer of ^1 to the large scales, thereby lending strong 
support to the prediction of no cascades. The large-scale energy spectra are 
shallower than k^^, and the small-scale spectra are steeper than k~^, consistent 
with the theoretical prediction. 
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